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We study analytically the correlations between the positions of tagged particles in the random 
average process, an interacting particle system in one dimension. We show that in the steady state 
the mean squared auto-fluctuation of a tracer particle grows subdiffusively ag (t) ~ t 1//2 for large time 
t in the absence of external bias but grows diffusively <Xq (i) ~ t in the presence of a nonzero bias. 
The prefactors of the subdiffusive and diffusive growths as well as the universal scaling function 
describing the crossover between them are computed exactly. We also compute cr%(i), the mean 
squared fluctuation in the position difference of two tagged particles separated by a fixed tag shift r 
in the steady state and show that the external bias has a dramatic effect in the time dependence of 
a'r(t). For fixed r, o1(t) increases monotonically with t in absence of bias but has a non-monotonic 
dependence on t in presence of bias. Similarities and differences with the simple exclusion process 
are also discussed. 



PACS numbers: 64.60.-i, 05.70.Ln 
I. INTRODUCTION 

Interacting particle systems in one dimension are 
amongst the simplest examples of many body systems 
that are far from equilibrium One of the most stud- 
ied examples is the simple exclusion process in one di- 
mension. In this system, each site of a one dimensional 
lattice is either occupied by a hardcore particle or it is 
empty. In a small time interval dt, each particle attempts 
to hop to the neighboring lattice site on the right with 
probability pdt, to the left neighboring site with proba- 
bility qdt and stays at the original site with probability 
1 — (p + q)dt. An attempted hop is completed provided 
the target site is empty. A wealth of results are known 
for this system Jl]-[|]. 

Another interacting particle system in one dimension 
that has attracted recent interest is the random average 
process (RAP) H|). In the RAP, particles are located on 
a real line as opposed to a lattice in the simple exclusion 
process. Let Xi(t) be the position of the i-th particle at 
time t (see Fig. 1). 
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1. The stochastic moves in the RAP. 



In a small time interval dt, each particle jumps to the 
right with probability pdt by an amount rf(xi+i — Xi), 
to the left with probability qdt by an amount r~ (xi — 
and stays at its original location with probabil- 
ity 1 — (p + q) dt. Here rf and r~ are independent ran- 
dom variables drawn from the interval [0, 1] with identical 
probability density function (pdf) f(r). Thus the jumps 



in either direction is a random fraction of the gap to the 
nearest particle in that direction. For convenience, we 
have defined the RAP with random sequential dynam- 
ics, though it has been studied with parallel dynamics as 
well [|||]. The detailed study of the RAP is important 
since it has shown up either directly or in disguise in a 
variety of problems including traffic models Q, models 
of mass transport || , models of force fluctuation in bead 
packs H , models of voting systems |?],[| , models of wealth 
distribution and the generalized Hammersley process 
p0| . Like the simple exclusion process, some aspects of 
the RAP are analytically tractable pP,pl|. In this paper, 
we derive some new exact results on the tracer fluctua- 
tions in the RAP where the dynamics of tagged particles 
are followed. 

The tracer diffusion has been studied in detail for the 
simple exclusion process and many interesting results are 
known [Q. In the exclusion process, the combined ef- 
fect of hardcore interaction and the external bias (p — q) 
shows up rather dramatically in the asymptotic long 
time behavior of the mean squared auto-fluctuation in 
the position of a tracer particle in the steady state. If 
Ci(t) = Xi(t) — (xi(t)) denotes the deviation in the po- 
sition Xi(t) of the i-th particle from its average value, 
then the mean squared auto-fluctuation is defined as, 
al{ta,U) + t) = ([&(to+t)-(i(t )] 2 ), where t is the wait- 
ing time after which one starts measuring the fluctua- 
tions. In the steady state to — ► oo, the asymptotic behav- 
ior of er 2 (t) = limto—i-oo cr 2 (io, + t) for large t is known 
Ej. In absence of external bias (p = q = 1/2), i.e. for 
the symmetric exclusion process (SEP), o"o(t) ~ At 1 / 2 for 
large t where the constant A = {2/-k) 1 / 2 (\ — p)/p is known 
exactly in terms of the density p of the particles [|l2] [l4| . 
This slow subdiffusive growth is due to the caging effect 
arising from hard core exclusion in one dimension where 
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a particle is always hemmed in by its neighbors. How- 
ever in the asymmetric case (ASEP) when a nonzero bias 
p — q > is switched on, one finds, somewhat unexpect- 
edly, ctq (t) ~ Dt for large t where the tracer diffusion 
coefficient D = (p — q)(l — p) fll5||l(f ]. The crossover 
from the subdiffusive to diffusive behavior of <r 2 (t), as 
an infinitesimal bias is switched on, was understood in a 
physically transparent way via a rather unusual mapping 
of the exclusion process to a (1 + l)-dimensional inter- 
face model |i~7 18 1. This mapping also established that an 
appropriately defined sliding tagged-particle correlation 
function varies anomalously as t 2 ^ 3 ]l7fl . This anomalous 
t 2 / 3 growth also shows up in the mean square fluctuation 
of the center of mass of the particles when viewed from 
a special moving frame flS[| , 

A question then arises naturally: what are the cor- 
responding results on the tracer diffusion for the RAP? 
The only known result is for the fully asymmetric RAP 
with q = (and time rescaled by p) where the particles 
move only to the right. In this limit, cr 2 (t) was computed 
by Krug and Garcia using a phenomenological hydrody- 
namic Langevin equation based on heuristic arguments 
as well as using an independent jump approximation 
Their result shows that a 2 (t) ~ D\t for large t with 
Di = p~ 2 — \i2) where p is the density of the par- 
ticles and pk = Jq drr k f(r) is the fc-th moment of the pdf 
/(r). Later, Schutz attempted to derive this result rigor- 
ously |20) by writing down the exact equation of evolution 
of the equal time correlation function G r (t) — (Co(OCr(0) 
and then using a chain of arguments. Note that the defi- 
nition <r 2 (t ,t + t)= ([0(to +t) — d(to)] 2 ) involves both 
the variance (( 2 (t)) which is an equal time observable 
as well as the unequal time correlation (d(to)(i(to + t)). 
Thus a proper approach, as followed in this paper, would 
be to compute these correlation functions exactly and 
then take the steady state to — > oo limit. 

The main results of this paper can be summarized as 
follows: 

1. We compute exactly the mean squared auto- 
fluctuation in the displacement of a single tracer 
particle, a 2 (t ,t Q + t) = ([Ci(*o + *) - C^o)] 2 ) for 
large to and t for all values of p and q in the 
RAP. In the steady state to — > oo, we show that 
al(t) = ]im to ^oc(7o{to,to + *) ~ A S R A pt 1/2 for 
large t for the symmetric RAP (SRAP) with p = q. 
For the asymmetric RAP (ARAP) where p > q, 
we find <r 2 (t) ~ Darap^ for large t. The con- 
stants A S rap = 2p- 2 (piJ, 1 /ir) 1 / 2 p, 2 /(pi - P2) and 
Darap = P~ 2 (p- q)ViP2/{pi - P2) are computed 
exactly. For the special case q = and p = 1, 
Darap reduces to D\ computed earlier in Refs. 



from the subdiffusive t 1 / 2 growth to the diffusive 
t growth as one switches on an infinitesimal bias 

(p-q)- 

3. We generalize the single tracer particle fluctuation 
(Tq (to , to + t) to the fluctuation in the position dif- 
ference of two tagged particles defined as <J 2 (t , t + 
t) = (Ki+r{to + t)- Ci(to)] 2 )- We show that in the 
steady state a 2 (t) = lim^^oo <r 2 (t n , t + t) grows 
monotonically with t for a fixed tag shift r for the 
SRAP. For the ARAP on the other hand, it grows 
with t in a non-monotonic fashion with a single 
minimum at a characteristic time t* = r / p\{p — q). 

4. We also compute various scaling functions that de- 
scribe the crossover of the tracer fluctuations from 
their non-steady state behavior to the steady state 
behavior as the waiting time to — > 00. 

The paper is organized as follows. In Sec. II, we de- 
fine the model precisely and set up our notations. In 
Sec. Ill, we calculate the equal time correlation function 
for the RAP for all p and q. Sec. IV contains the ex- 
act calculation of the unequal time correlation function. 
In Sec. V we compute the mean squared fluctuation in 
the displacement of a single tracer particle. The Sees. 
VA and VB contain respectively the discussions on the 
SRAP and the ARAP, while the crossover between them 
is discussed in Sec. VC. The Sec. VI contains the gener- 
alization to the two-tag correlation functions. Finally we 
conclude with a summary and discussion in Sec. VII. 



II. THE MODEL AND PRELIMINARIES 

We consider a system of particles of average density p 
located on a real line. Let Xi(t) denote the position of 
the i-th particle at time t (see Fig. 1). In an infinites- 
imal time interval dt, each particle jumps with proba- 
bility pdt to the right, with probability qdt to the left 
and with probability 1 — (p + q)dt it rests at its original 
location. The actual amount by which a particle jumps 
(either to the right or to the left) is a random fraction 
of the gap between the particle and its neighboring par- 
ticle (to the right or to the left). For example, the jump 



to the right is by an amount r^(xi + i 



and to the 



left by r~(xi — Xi-\). The random variables rf are in- 
dependently drawn from the interval [0, 1] and each is 
distributed according to the same pdf f(r) which is arbi- 
trary. We start from an arbitrary but fixed initial condi- 
tion at t = and averaging of physical quantities is done 
over all histories of evolution keeping the initial condition 
fixed. The time evolution of the positions a;j(t)'s can be 
represented by the exact Langevin equation 



2. We compute exactly the universal scaling func- 
tion that describes the crossover behavior of a 2 (t) 



Xi(t + dt) = Xi(t) + ji(t), 
where Ji(t) are random variables given by 



(1) 



2 



!r^~(x,+i(t) — Xi(t)) with prob pdt, 
r~(x l ^ 1 (t) - Xi(t)) with prob qdt, (2) 
with prob 1 — (p + q)dt. 

The random variables rf 1 are independent and each is 
distributed over the interval [0, 1] with the same pdf 
/(r). The k-th moment of the pdf is denoted by fik = 
J a drr k f(r). Note that since < r < 1 and f(r) > 0, 
Pi > P2- 

We define a new random variable which measures 
the deviation of Xi(t) from its mean value 



Ci(t) = Xi{t) - (Xi(t)). 



(3) 



From Eqs. (|l|) and (|^), one can easily derive the evolution 
rules for the Q variables. We find 

&(t + <ft) = Ci(*) - (p - ?)— * + ( 4 ) 



where ^i(i) is given by 



»h(t) = 



T(Ci- 



■i(*) - d(t) + p- 1 ) 
-i(t) - d(t) - p- 1 ) 



o 



with prob pdt, 
with prob qdt, (5) 
with prob 1 — (p + <?)dt. 



By definition, (&(*)) = 0. Also from Eq. (g), it follows 
that =(p-q)fdt. 

In this paper, we will focus on the mean squared dis- 
placement of a tagged particle. It turns out that the 
asymptotic behavior of the mean squared displacement 
depends crucially on whether one starts measuring these 
fluctuations after some finite waiting time to or if one 
first waits for an infinite time and then starts measur- 
ing the statistics. The latter corresponds to measuring 
the fluctuations in the steady state. This is similar to the 
'approach to stationary' versus 'stationary' regimes found 
in various interface models This can be quantified 
precisely in terms of the following correlation function, 

<rt(to,t Q +t) = ({ti(t + t )-ti(t )) 2 ), (6) 
= G (t+t ) + Go(to) ~ 2C (t ,t +t), (7) 

where G r {t) = (d(t)(i +r (t)) is the equal time correla- 
tion function and C r (to,io + 1) = (d(to)Ci+r(to + 1) with 
t > denotes the unequal time correlation function. For 
t = 0, the unequal time correlation function reduces to 
the equal time correlation function, C r (to,to) — G r (to). 
Note that we have assumed an infinite system so that the 
translational invariance holds. In the next two sections 
we calculate analytically the correlation functions G r (t) 
and C r (to, to + t) respectively. 



III. EQUAL TIME CORRELATION FUNCTION 

In this section, we calculate the equal time correlation 
function G r (t) = (Q(t)Q +r (t)) exactly for the RAP for 



all p and q. Our starting point is Eq. (Eh in conjunction 
with Eq. (||) describing the evolution of the £i variables 
with time. We consider the evolution equations (Eq. (H)) 
for both Q(t + dt) and Ci+r(t + dt), multiply them and 
then take the average () over all histories, keeping terms 
only upto 0(dt). This yields, in the limit dt — > 0, the ex- 
act evolution equation of the correlation function G r (t) 
and we obtain, 



jG r {t) = m{p + q) [G r+1 (t) + G r _i(t) - 2G r (t)} 

+ <W 2 (P + <Z) [ P - 2 + 2{G (t)-G 1 (t))] . (8) 

The Eq. (H) is valid for all positive and negative inte- 
gers r including r = and clearly G r (t) = G- r (t). Thus 
the equation of evolution for the two point correlations 
involve only two point correlations and not higher order 
correlations. This closure property is crucial for obtain- 
ing an exact solution for the correlation functions. The 
key reason behind this closure lies in the fact that the 
random fractions r^'s at time t are independent of the 
Q(t). One noteworthy fact about Eq. (||) is that the 
rates p and q make their appearance only as an overall 
multiplicative factor (p + q). We could absorb this factor 
into the time by doing a suitable rescaling, and hence, the 
equal time correlation function G r (t) is same for both the 
ARAP and the SRAP. 

We note that this equation was also derived in Ref. 
pot by a rather lengthy method, but was left unsolved. 
In this section, we derive an exact solution of Eq. (^). 
Note that even though Eq. (||) represents the diffusion 
equation (in discrete space) with a source term at the 
origin r = 0, its solution is nontrivial due to the fact 
that the source term depends on Go(t) and G\(t) which 
need to be determined self-consistently. Similar diffusion 
equations with source term for the correlation functions 
have also appeared recently in the context of aggregation 
models with injection |2^| . Before proceeding to solve 
Eq. (||), we first set up our notations. We define the 
standard Fourier transform 



G(k,t) = G r(*)e 



ikr 



the Laplace transform 



G r (s) = / G r (t)e- St dt, 
Jo 



(9) 



(10) 



and the joint Fourier-Laplace transform, 

POO °° 

F{k,s)= / G{k,t)e- St dt= V G r {s)e lkr . (11) 
Jo „_ „ 



r— — oo 



Taking the joint Fourier-Laplace transform of Eq. (j^) 
we obtain 
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F(k,s) 



m{p + q) [ P - 2 + 2s (Go(a)-Gi(a)) 



s[s + 2/ii(p + - cos A:)] 



(12) 



where we have assumed that initially G r (0) = which 
is true for any fixed initial condition. For random initial 
condition, F(k, s) will contain additional terms arising 
from the initial condition, but one can show that they do 
not contribute to the asymptotic large time properties 
of G r (t) as long as the initial condition has only short 
ranged correlations. We therefore use G,(0) = without 
any loss of generality 

The Eq. (|l2|) contains two unknowns Go(s) and Gi(s). 
One of them, say Gi(s) can however be expressed in 
terms of Gq(s) by taking directly the Laplace transform 
of Eq. (|) for r = and using Gi(t) = G_i(f). This 
gives the relation 



2(/xi-/i 2 ) (<3o(a)-Gi(a) 



sGq(s) = (p+q) 



Substituting Eq. @ in Eq. (|T|) we obtain 

pi{p + q)p~ 2 - s 2 G (s) 



F(fc,s) 



/'2 



(Mi - M2) s [s + 2/ii(p + g)(l - cosfc)] ' 



(13) 



(14) 



We now have to determine Gq(s) self-consistently. This 
can be done by using the inverse Fourier transform 



G r (s) 



1 

2^ 



F{k,s)e- lkr dk 



(15) 



Substituting the expression of F(k, s) from Eq. ( O ) in 
Eq. ( |l5| ) at r = 0, we obtain the exact Gq(s) 



G (s) = 



/ii/i 2 (p + g) 



p- 2 /(0, S ) 



(Mi - M2) 



+ / s 5/(0,s) 
0*1 -/*a) v ' ; 



(16) 



where 7(r, s) is given by the integral 

e~ lkr dk 



'<••••> "s£ 



[s + 2// x (l - cosfc)] 



'2f^ +S- y/fi 2 + 4gj 



(17) 



where Mi = Mi(p + <7)- Knowing Go(s) determines F(fc, s) 
completely by Eq. (|T4|) and hence G r (s) for all r by the 
Fourier inversion formula in Eq. (O). We obtain 



G r (s) 



MiM2(p + ?) 



i 2 I(r,s) 



(Mi - M2) 



where 7(r, s) is given by Eq. (fi"7|). 



s/(0,s) 



(18) 



To obtain G r (t) we need to perform the inverse Laplace 
transform G r {t) = £~ 1 {G r (s)] with respect to s. In gen- 
eral for arbitrary t this is difficult. However, for large t, 
this inverse can be obtained in closed form. For large t, 
one needs to consider the small s behavior of G r (s) in 
Eq. (fi8|). Let us first consider the case r = 0. Putting 
r = in Eq. (|l7|) and taking the s — > limit we find to 
leading order, 



J(0,s) 



1 



(19) 



2v / Mi(P + <?)s 

Substituting this small s expression of 1(0, s) in Eq. © 
and taking the inverse Laplace transform we find that to 
leading order for large t, 



G (t) 



y /fJ, 1 {p + q)fJ, 2 p 
(Mi - M2)V7r 



-Vt. 



(20) 



Next we consider the behavior of G r (t) for |r| > 0. 
From Eq. (|l^) it is clear that the appropriate scaling 
limit consists of taking the limit s — > 0, \r\ — ► 00 but 
keeping |r|i/s fixed. In this scaling limit, Eq. ( |l7j ) yields, 



Z(r,s) 



f -H^ 

2v / Mi(P + '?) s ^ I \/Mi(P + 9) 



exp 



(21) 



We note that the formula for I(r, s) in Eq. ( pl| ) reduces 
to Eq. ( p^ ) for |r| = 0. This indicates that even though 
Eq. (2l|) was derived in the scaling limit, it continues to 
hold even for r = 0. 

Substituting this small s expression of I(r, s) in Eq. 
( |l8| ) and taking the inverse Laplace transform we obtain 
for large t, 



G r (t) 



\Jpi{p + q)p2p" 
2(mi - M2) 



-£- 



-3/2 p-M^/s/ii^ip+q)] 



(22) 



Fortunately the inverse Laplace transform in Eq. ( |2z 
can be done in closed form which gives us the following 
asymptotic scaling behavior of the equal time correlation 
function G r (t), 



(Mi - M2)V^ 



(23) 



Here f± (y) is a universal scaling function independent of 
the model parameters such as p, q and the moments Mfc 
of the pdf f(r) and is given by 



v - VtF y erfc(y), (24) 

— 2 

e~ u du is the standard com- 



where erfc(y) = 
plimcntary error function. This scaling function has the 
asymptotic behaviors, fi{y) *~ 1 — \f^y as y — > and 
~ y~ 2 e~ v /2 for y — > 00. 
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As a final remark, we note again that if one puts \r\ = 
in the formula for G r {t) in Eq. (^J) one recovers the cor- 
rect Go(t) as given by Eq. @. Thus the scaling range 
includes even the r = point. The Eq. ( |20| ) thus pro- 
vides us the exact behavior of the first two terms in the 
expression for <7g(to,to + t) m Eq. (Q). The remaining 
task is to evaluate the third term in Eq. (Q) which in- 
volves the unequal time correlation function and this is 
done in the next section. 



IV. UNEQUAL TIME CORRELATIONS 

In this section we compute the two time tag-tag cor- 
relation function C r (t 0: t + t) — (d(to)Ci+r(to + t)) for 
the RAP. We start at time to and then evolve the Ci+r 
variables by Eq. (0) for all subsequent time. Let us first 
rewrite the Eq. (0) at time to + t + dt, 

(i+r{t + t+dt) = ( i+r (t + t) - (p-q) — dt +7] l+r (t + t). 

P 

(25) 

We then multiply both sides of Eq. ( p5| ) by Ci(^o) an d 
average over the noise keeping terms only upto 0(dt). In 
the limit dt — > 0, we obtain the exact evolution equation 
of the two time correlation function, 



dCr(t ,t +t) 
dt 



Hi [pC r +i(t Q ,t Q + t) + qC r -i(t Q ,t + t) 
- (p + q)C r (t , t + t)} for t > 0. (26) 



Note that at t = 0, the unequal time correlation 
function reduces to the equal time correlation function 
C r (to,to) = G r (to). Thus starting at t = with 
the initial condition C r (to,to) — G,-(to), the function 
C r (to,to + t) evolves with time t according to the Eq. 

As in the preceding section, we define the Fourier 
transform C(k,t ,t + t) = E^-™Cr(*o,*o + t)e ikr . 
Taking the Fourier transform of Eq. (pfl) we obtain 



C(k,t ,t +t) = G(Mo)e- WQ(fe)t , 



(27) 



where a(k) — p + q — (pe 



qe 



ik) 



and G(k, to) is the 



Fourier transform of the equal time correlation function 
as defined by Eq. (H) . Taking further the Laplace trans- 
form H(k, s, t) = ff C(k, t ,t + t)e- st °dt of Eq. (§|) 
we obtain 



H(k,s,t) = F(k,s)e- f " ia ^\ 



(28) 



where F(k,s) is given exactly by Eq. (14) with Gq(s) 
determined from Eq. ([l6|). 

Proceeding as in the previous section, the Laplace 
transform C r (s,t) — J °°C r (to,io + t)e~ st °dto can then 
be determined from the joint Fourier-Laplace transform 
H(k, s,t) by the inversion formula 



C r (s,t) = — 

Z7T 



H(k, s, t)e 



dk, 



(29) 



where H (k, s, t) is given by Eq. (p8|). Substituting in Eq. 
( ^9| ) the exact expression of F(k, s) from Eq. ((IJ) and 
that of Go(s) from Eq. (|l6|), we obtain the following final 
expression of the Laplace transform 



C r (s,t) 



piMp + q) 



(Mi ~ Ma) 



1 



1 



e —ikT—ii-i_a.{h)t^u 



sl{0,s) 



2tt [s + 2m (p + q) (1 - cos A;)] ' 



(30) 



Note that for t = 0, CV(s, t) as given by Eq. ( J30| ) reduces 
to G r (s) given by Eq. ( |lq) as expected. The equation 
( ^0| ) is central to our subsequent analysis for various lim- 
iting behaviors. 



MEAN SQUARED TRACER 
AUTO-FLUCTUATION 



In this section we calculate <jQ(to,to + t) in the RAP 
using the exact results for the equal time and two time 
correlation functions obtained in the previous sections. 
We consider first the symmetric case SRAP with p = q 
in subsection A followed by the derivation for the asym- 
metric case ARAP with p > q in subsection B. In sub- 
section C, we show how the steady state fluctuation 
(Tq(£) = limta—KxCTafto, to + t) crosses over from the subd- 
iffusive behavior to the diffusive behavior as one switches 
on an infinitesimal bias and we calculate the crossover 
scaling function exactly. 



A. SRAP 

Here we consider the symmetric case p = q. For the 
calculation of CTq(£o, *o + t) we only need the asymptotic 
behavior of C r (to,to + t) for r = as evident from Eq. 
(0). To obtain Cb(io> to + t) we need to invert the Laplace 
transform in Eq. ( |30"| ) for r = and p = q. As before, 
this inversion is difficult in general for all to. However 
the finite but large to limit can be worked out by analyz- 
ing the small s behavior of Eq. (|30|). It turns out that 
the appropriate scaling limit in this case involves taking 
s — > 0, t — > 00 but keeping st fixed. In this scaling limit, 
the integration in Eq. ( |30| ) can be carried out in closed 
form and we obtain (with p = q), 



C (s,t) 



2(/xi - /x 2 )s 3 / 2 



(31) 



We then need to invert the Laplace transform in Eq. ( |3l| ) 
with respect to s to obtain the asymptotic behavior of 
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Co(to,to + t) for large to- Fortunately this inversion can 
be done in closed form and we obtain 



Co(to, to + t) 



(mi - m)V^ 



Vhf2 



t 

2k> 



(32) 



where the scaling function f 2 (y) is again universal and is 
given by, 



h(y) = V T +y-Vv- 



(33) 



We are now ready to compute erg (to, to + t) from Eq 
( 71) . Using the result for the equal time correlation in Eq 
(2w and the one for the two time correlation in Eq. (32) 
we obtain from Eq. (Q) our main result 



o"o(*0) to + t) — 
(Ml - A*2)Vtt 



t 

2T Q 



(34) 



where fq(y) is given by Eq. fl33j). Note that this result 
in Eq. (Q) is derived in the scaling limit when both to 
and t are large with their ratio t/to kept fixed. 

We now discuss two different limits of Eq. ([34]). First 
we consider the steady state limit to — ► oo with t large 
but fixed. In this limit, Eq. (E34J) yields 



o-o(t) = limto^ooO'o(*o,*o + t) 



(Mi ~ M2)v^ 



\/t. (35) 



In the opposite limit, when the waiting time to is finite 
(away from the steady state) but the evolved time t goes 
to infinity, we obtain from Eq. (|34|, 



Iim 4 ^ 00 o'o(to,io + t) 



y/2pp 1 p 2 P 

(Ml - M2)V^ 



(36) 



Thus the mean squared auto-fluctuation in these two op- 
posinglimits differ by a factor y/2. The Eqs. (g), @ 
and ( p6|) are amongst the important new results of this 
paper. 



B. ARAP 

In this subsection we calculate 0o(to,to + t) in the 
asymmetric case when p > q. Once again we have to 
invert the Laplace transform in Eq. (^0|) for r = but 
now with p > q. In this case it turns out the appropriate 
scaling limit consists of taking s — > 0, t — > oo as in the 
SRAP but keeping -y/it instead of the scaling variable si 
in the SRAP. In this scaling limit, the integration in Eq. 
© with r = yields 



C o(M ) = M +y e -(HV^W* (37) 
2(mi - P2)s i ' 2 

The Laplace transform in Eq. (|3^) can be inverted as in 
Eq. (E2|) and we obtain 



Co(*o, *o + 1) — — fpS^ — 1^ 2 P — v ^tof 1 

(Mi - M2)V7T 



(38) 



where the universal scaling function /i(y) = — 
y/iry erfc(y) is the same as in Eq. (p3). 

Substituting the results in Eq. ( B8| ) and Eq. (|20| ) in 
Eq. (0) we obtain 



2/. , i i\ y / Mi(P + g)M2P 
(Mi - M2 V 71 " 



m(p-q)t 



2^{p + q)t 



(39) 



As in the SRAP we now discuss the two different lim- 
its. In the steady state t — > oo with fixed large t we 
obtain from Eq. (p9|), 



<7q(*) = hm to _ >00 a-o(to,to + *) = 



M1M2M 2 (p-q) 
(mi - M2) 



t. (40) 



Thus in this case a^t) grows diffusively for large t, 
Uo(t) = D ARApt where the diffusion constant, 



Darap = P 2 {p- q) 



M1M2 
(Mi ~M2) : 



(41) 



depends explicitly on p and q. For q = and p = 1, it 
reduces to the expression D\ — p~ 2 M1M2/ '(Mi — M2) de- 
rived by Krug and Garcia using the independent jump 
approximation jj] and later rederived by Schiitz |20| us- 
ing a different approach. 

We make a brief comment here on the approach used in 
Ref. El in deriving the diffusion constant D\. In his ap- 
proach, Schiitz started with the evolution equation (||) for 
the equal time correlation function and then used a chain 
of arguments to derive the diffusion constant D\. His ap- 
proach didn't require any knowledge of the two time cor- 
relation function or even the solution of the equal time 
correlation function. As evident from the definition in 
Eq. (0) that <Tq (to ,to+t) requires the knowledge of both 
the equal and the two time correlation functions. Thus 
it was rather remarkable that the correct value of the 
diffusion constant for q = and p = 1 was recovered in 
Ref. pQ]. However this turns out to be purely fortuitous. 
Note that the evolution equation (||) is independent of the 
bias in the system. Thus the approach of Schiitz would 
predict that the diffusion constant is also completely in- 
dependent of the bias (p — q) and is always given by D\ 
(provided t is scaled by (p + q)). This is clearly wrong 
as evident from the exact expression in Eq. (pd|). In 
particular for the symmetric case p — q = 1/2, the ar- 
guments of Ref. pi would predict a diffusive growth of 



G 



cTq (t) with the diffusion constant D\ . This is again incor- 
rect since for p = q the diffusion constant is from Eq. 
( fi"l| ) which is consistent with the correct asymptotic sub- 
diffusive growth of Oo(t) as given exactly by Eq. (35). 
The problem in the derivation of Schiitz can be traced 
back to the fact that his arguments only used equal time 
correlations (which involve only (p + q)) and not the two 
time correlations. The dependence on the bias (p — q) of 
the diffusion constant Darap comes only from the two 
time correlations. The derivation of Ref. |2(]] misses this 
important fact and is rather fortuitous to obtain the cor- 
rect value Di of the diffusion constant for the special case 
when p = 1 and q = 0. 

We end this subsection by discussing the other limit 
when the system is away from the steady state, i.e. when 
to is large but finite and t — > oo. In this limit, we obtain 
fromEq. (f§) 



limt^ooCToCto.io +t) 



V 'pi{P + q)p2R 
(Mi ~ m)V^ 



-Vt, (42) 



the same result as in the SRAP in this limit [Eq. (p6[)]. 
Thus away from the steady state the tracer particle 
doesn't sense the presence of bias. The exact result in 
Eq. ( f42| ) is consistent with that of Krug and Garcia us- 
ing a phenomenological hydrodynamic equation 8 . 



C. Crossover Between SRAP and ARAP 

In the previous subsections, we have seen that the 
asymptotic large t behavior of a 2 (to, to + t) does not 
depend on the bias (p — q), when the system is away 
from the steady state (finite to)- However, in the steady 
state (to — > oo) it behaves rather differently in the sym- 
metric and asymmetric cases. In the steady state of the 
SRAP (p = q), crg(t) - t 1 / 2 while for the ARAP (p>q), 
Oo(i) ~ t. Thus a natural question is: How does the be- 
havior of <7q (t) crosses over from the subdiffusive growth 
for p — q to the diffusive growth as one switches on an 
infinitesimal bias (p — q)l In this subsection we compute 
exactly the universal scaling function that describes this 
crossover behavior of <t 2 if) ■ 

To calculate the crossover behavior we return to our 
central equation ( |30| ) with r = 0. We have seen in the 
previous subsections that in the scaling limit s — > and 
t — ► oo of the Eq. (|3(]) , the appropriate scaling variable 
that is kept fixed is st for p = q, where as, it is i/st for 
p > q. Thus, to compute the crossover behavior, we need 
to keep the leading order terms in both of these scaling 



variables fixed while expanding Eq. (30) for small s and 
large t. This makes the calculation of the crossover be- 
havior somewhat delicate. To leading order, we find after 
elementary algebra 



C (s,t) 



y / Pl(p + q)p 2 p 
(pi - M2)s 3 / 2 



2tt 



-i{p-q)y/ His/(p+q) tz-stz 2 /2 



1 + Z 2 



-dz. 



(43) 



Note that for the symmetric case p = q, the integral in 
Eq. ( |4g| ) can be done and we get back Eq. ( j3l|) of Sec. 
VA. Similarly, for the asymmetric case p > q, in the limit 
s — > keeping the scaling variable y/st fixed, one drops 
the second term in the exponential in the integrand of Eq. 
d43| ) and performing the resulting integral we recover the 
Eq. @ of Sec. VB. 

To compute the crossover behavior we need to keep 
both the terms inside the exponential in the integrand of 
Eq. ( ff3| ) and perform the integral. Fortunately this inte- 
gral can be done in closed form using the standard con- 
volution theorem. We omit the details here and present 
only the final result, 



C (s,t) 



y/pi(p + q)p 2 p" 

4(mi - M2)s 3/2 
(2u- 



' J erfc 



V 2^ 



erfc 



2u 



2Vu 



(44) 



where u = st/2 and v — (p — q)\/pis/ [p + q) t. We then 
expand the Eq. (|44|) further for small s to obtain the 
steady state to — > oo behavior. Note that we needed 
to first do the integral in Eq. (^3|) and then take the 
s — > limit. The reverse order unfortunately does not 
work. Expanding Eq. (Q) for small s, keeping only the 
leading order terms in s and finally inverting the Laplace 
transform of the resulting expression we obtain for large 
to 



Co(toi to + 1) 



{pi - 



P2)V 7! ' 

^/¥pl(p - q)t 



Z^pTq 



erf (w(t)) 



(45) 



where w(t) = (p- q)y / pit/ [2(p + q)]. 

We now use the results from Eq. (|45|) and Eq. (^0|) in 
Eq. (^) and eventually take the strict to — > oo limit to 
obtain the final form of the steady state auto-fluctuation 



al(t) = limt o -KX)0o(*o,£o +t) 



pip 2 (p-g)p~ 
(Pi ~ Pi) 



-tY 



V 2 (p + v) 



(46) 



where Y(y) is a universal crossover scaling function given 
by 



Y(y) =er%) + 



1 e~ v 

fx y 



(47) 



The scaling function has the asymptotic behavior Y(y) ~ 
l/( v / 7rj/) as y — * and Y(y) — > 1 and y — * oo. Note that 
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for fixed p — q > 0, if we take the limit t — > oo in Eq. ( fig ) 
(which corresponds to y — > oo in the scaling function in 
Eq. ([47j)) we recover the result of Eq. Similarly 
if we take the p — q — > limit for fixed i in Eq. 
(corresponding to taking y — * limit in the scaling func- 
tion Y(y)), we recover, as expected, the result of Eq. 
( p5| ) of the symmetric case. Thus the Eq. (ff6| ) and the 
associated scaling function F(y) in Eq. ([47[) describes 
the crossover behavior from the subdiffusive to diffusive 
growth as one switches on an infinitesimal bias. 



VI. GENERALIZATION TO THE TWO- TAG 
CORRELATION FUNCTION 

So far in this paper we have concentrated only on 
the mean squared auto-fluctuation of a tracer particle, 
Oq(*o,*o +t) = {(d(t + t ) - G(*o)) 2 )- A natural gener- 
alization of the auto-fluctuation would be to study the 
two-tag correlation function defined as 

a?,(t ,t + t) = ((Q+ r (t + to)-Q(to)) 2 ), (48) 
= G (t+t ) + Go(io) ~ 2C r (t ,t +t), (49) 

where G r (t) and C r (to,to + t) are the usual equal time 
time and the two time correlation functions already de- 
fined and derived in the previous sections. Note that for 
r = 0, the two-tag correlation in Eq. (^8|) reduces to the 
single tag function aft (to, t + t). 

Of particular interest would be to compute the two- 
tag correlation function in the steady state, i.e. 07 (i) = 
limt _KX3 cjKto, to+t). For the exclusion process this two- 
tag correlation function was first introduced in Ref. fl§l l 
and the presence of bias was found to have a dramatic 
effect on the time dependence of af(t) for a fixed r. It 
was found numerically that while in the SEP crj?(t) in- 
creases monotonically with t for a fixed tag-shift r, in 
the ASEP <r;?(i) has a non-monotonic dependence on t 
fill . In the ASEP of (i) first decreases with time t, be- 
comes a minimum at some characteristic time t* and then 
starts increasing again. A harmonic model was proposed 
in Ref. [is] for which af(t) could be computed analyti- 
cally and was found to be in qualitative agreement with 
the numerical results of the exclusion process. But to the 
best of our knowledge, exact calculation of <r^(t) for the 
exclusion process is still an unsolved problem. However 
it turns out that for the RAP, it is possible to compute 
this function of(i) exactly for large t. The exact solution 
of a^(t) in the RAP, as shown below, shares the similar 
features as in the exclusion process. 

From Eq. (^9|) it is evident that we just need to com- 
pute the large to behavior of the two time correlation 
function C r (to, to + t) for fixed nonzero r. In the previ- 
ous sections we have analyzed in detail the r = case. 
It turns out that the analysis for r ^ proceeds more or 
less in the same manner as in the r = case. We start, 



once again, from the central equation (|3(j) . To avoid sep- 
arate calculations for the SRAP and the ARAP, we take 
the line of approach used to calculate the crossover be- 
havior in subsection VC. For r 7^ 0, it turns out that the 
equation ([43|) gets replaced by a similar looking equation, 



C r (s,t) 



"2tt 



yVl(p + q)H2P 

(a*i - A«2)s 3 / 2 



00 e -izy/s/[m(p+q)]R-stz 2 /2 



1 + Z 2 



-dz, 



(50) 



where R = r + /zi {p — q)t signifies the drift of the parti- 
cles to the right with average velocity ni(p—g) for p > q. 
Clearly for r = 0, Eq. ( |50| ) reduces to Eq. ((43[). Starting 
with Eq. (^(j|) we then follow exactly the same steps as 
used in subsection VC. Since the steps are identical we 
skip all the details and present only the final result. In 
the strict steady state limit t — > 00, we finally obtain 
the following scaling form 



al(t) = lim ta -,oocrr(to,to + t) 



y / 2fi 1 (p+q)Li 2 p~ 
(Ml - A*2)V^ 



-VtW{ 



R 



^2 ( l 1 {p + q)t 



(51) 



where R = r + (p — q)t and W(y) is again a universal 
scaling function given by, 



W(y) 



~ y2 + \/7n/erf(y). 



(52) 



Clearly Ui(p — q)t represents the average drift while 
l{t) = yj2[i\{p + q)t represents the diffusive length scale. 

We note that the scaling function W(y) is a symmetric 
function of y about y — with a minimum at y — 0. For 
the SRAP, p — q and hence R = r. Thus for a fixed r, 
it follows from Eq. ( |5l|) that ofW increases monotoni- 
cally with t. For the ARAP on the other hand, p > q 
and R = r + fii(p — q)t. If one fixes r to a negative 
value and increases t, the variable R remains negative 
till the characteristic time t = t* = r//ii(p — q), be- 
yond which it becomes positive. The scaling variable 
y = R/ y/2/j,i(j> + q)t behaves in the same way. Thus 
<J 2 r {t) in Eq. ( |5l| ) first decreases with time, becomes a 
minimum at t = t* = —r/(ii(p — q) and then starts in- 
creasing again. In Fig. (2) we plot the function <rj?(i) 
in Eq. @) for both the SRAP (with p = q = 1/2) and 
the ARAP (with p — 1 and q = 0) for the same value 
of r = —2 and choosing the parameter values [i\ = 1/2, 
/U2 = 1/4, p = 1. These features in the RAP, derived here 
exactly, are qualitatively similar to those in the exclusion 
process studied in Ref. pi. 



8 
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FIG. 2. The steady state two-tag correlation function cr^(t) 
in Eq. ( |5l| ) plotted as a function of t for fixed r = — 2 for pa- 
rameter values fii — 1/2, fi2 = 1/4 and p = 1. The solid line 
shows the monotonic growth of cr^(t) with t for the SRAP 
(p = q = 1/2) while the dashed line shows the non-monotonic 
growth for the ARAP (p = 1 and g = 0) . 



VII. CONCLUSIONS 



In this paper we have studied analytically the mean 
squared fluctuations in the diffusion of both a single 
tagged particle and two tagged particles in the ran- 
dom average process (RAP) for all values of the hop- 
ping rates p and q in one dimension. We have shown 
that in the steady state, the auto-fluctuation of a tagged 
particle grows subdiffusively as cr§(t) ~ AsRApt 1 ^ 2 for 
p = q and diffusively o"o(i) ~ DARApt for p > q where 
A S rap = 2p- 2 (pni/ir) 1 / 2 n 2 /(pi ~ M2) and D A rap = 
P 2 {p ~ 1)^1^2/(^1 — aO- These behaviors of <7q(£) 
are similar to those in the simple exclusion process ex- 
cept the prefactors A = (2/vr) 1 / 2 (l - p)/p and 
D = (p — q)(l — p) ]l5| , [l6t are different in the exclusion 
process. Besides the steady state mean squared two-tag 
fluctuation a 2 (t) in the RAP grows monotonically with t 
for p = q and non-monotonically for p > q 1 in much the 
same way as in the exclusion process. 

These findings raise the question whether or not the 
RAP is in the same universality class as the simple ex- 
clusion process in one dimension. Perhaps the RAP is 
just a coarse grained version of the exclusion process in 
one dimension? The answer to this question seems to 



be in the negative due to a very crucial difference be- 
tween the two processes. In the exclusion process for 
p > q, it is well known that there exists an anomalous 
£2/3 g row th hidden in the problem apart from the usual 
t 1 / 2 and t growth pj| , ^7| . This anomalous growth shows 
up either in the mean squared fluctuation of the center 
of mass of the particles when viewed from a special mov- 
ing frame or alternately in the two-tag correlation 
function a 2 (t) if one chooses the tag shift r to be sliding 
with time with a special velocity r = —p 2 (p — q) 17|. It 
turns out that the prefactor of this t 2 ^ 3 growth is propor- 
tional to oc d where j(p) is the current density in a 



hydrodynamical description |17|. For the exclusion pro- 
cess, j(p) = p(l — p) and hence the prefactor is nonzero. 
For the RAP on the other hand, j(p) — pi(p — q) and 
is independent of p. This is because j(p) = p(v) where 
the average velocity (v) = p,\ (p — q)/p as can be easily 
derived from Eqs. (|l|) and (Q). As a result, for the RAP, 
the anomalous t 2 / 3 growth is absent which puts it in a 
different universality class than the simple exclusion pro- 
cess. In this sense the RAP seems to be closer to the 
harmonic model studied in Ref. [[isfl . 

In this paper we have considered the RAP only in one 
dimension. An obvious generalization would be to higher 
dimensions. A natural way to generalize the model to 
higher dimensions would be as follows. One considers 
particles located in the continuous <i-dimensional space. 
In a small time interval dt, each particle makes a list of 
all its nearest neighbors in various directions in space, 
chooses one of them at random and jumps in the corre- 
sponding direction by a random fraction of the Euclidean 
distance to that neighbor. This is an isotropic version, a 
generalization of the SRAP. Similarly one can define an 
anisotropic version as well. To the best of our knowl- 
edge, the RAP has not been studied so far in higher 
dimensions. The question of tracer diffusion in higher 
dimensions, especially in two dimensions where one may 
expect a logarithmic correction, also remains completely 
open. 
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